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Abstract. We present a preliminary local thermodynamic equilibrium (LTE) abundance 
analysis of the template halo red giant HD 122563 based on a realistic, three-dimensional 
(3D), time-dependent, hydrodynamical model atmosphere of the very metal-poor star. We 
compare the results of the 3D analysis with the abundances derived by means of a standard 
LTE analysis based on a classical, ID, hydrostatic model atmosphere of the star. Due to the 
different upper photospheric temperature stratifications predicted by ID and 3D models, 
we find large, negative, 3D- ID LTE abundance differences for low-excitation OH and Fe i 
lines. We also find trends with lower excitation potential in the derived Fe LTE abundances 
from Fe i lines, in both the ID and 3D analyses. Such trends may be attributed to the 
neglected departures from LTE in the spectral line formation calculations. 

Key words. Convection - Hydrodynamics - Line: formation - Stars: abundances - Stars: 
atmospheres - Stars: late-type - Stars: individual: HD122563 



1. Introduction 

HD 122563 is one of the most widely studied 
halo stars and is generally considered as a tem- 
plate for met al-poor red giants in abundance 
analyses (e .g., Barbuv et al.l2003l:ICowan et all 
l2005t lAokietal.1120071) . Up to now, all spec- 
troscopic analyses of this star have relied on 
the use of classical, one-dimensional (ID), sta- 
tionary, hydrostatic model stellar atmospheres. 
In this contribution, we present some prelimi- 
nary results from the first abundance analysis 
of HD 122563 based on a three-dimensional, 



time-dependent, hydrodynamical model atmo- 
sphere of the halo giant. 

2. Methods 

We use fhe stagger- code 

dNordlund & Galsgaardl Il995lfl to carry 
out a 3D, radiative, hydrodynamical, simula- 
tion of convection at the surface of a red giant 
with approximately the same stellar param- 
eters as HD122563, that is with an effective 



1 http : //www . astro . ku . dk/~ kg/Papers/ 
MHD_code.ps. gz 
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Fig. 1. Gray shaded area: Temperature distribution 
as a function of continuum optical depth at A = 
5000 A in the 3D metal-poor red giant surface con- 
vection simulation used in the present analysis of 
HD122563; darker areas indicate values with higher 
probability. Solid line: Temperature stratification in 
the ID marcs model atmosphere constructed for the 
same stellar parameters. 



temperature r e ffx4600 K, a surface gravity 
logg =1.6 (c gs), and a scaled solar che mical 
composition ( Grevesse & SauvaH 1 19981) with 
[X/H]= -3 E f° r all metals. The equations 
for the conservation of mass, momentum, and 
energy are solved together with the radiative 
transfer equation on a discrete mesh with 
480x480x240 numerical resolution for a 
representative volume of stellar surface deep 
enough to cover about eleven pressure scale 
heights and large enough to incorporate about 
ten granules (-3700x3700x1 100 Mm 3 ). 

We adopt open boundaries at the top and 
at the bottom of the simulation domain, and 
periodic boundary conditions at the sides. We 
use a realistic equation-of-state (Mihalas et al. 
1988b and conti n uous and li ne opaci t ies by 
Gusta fsson et"ai1 (1 19751) and iKuruczl (1 19921 
1993b . To compute the radiative heating rates 
entering the energy conservation equation, 
we group the o pacities in four opacity bins 
(lNordlundl ll982) and solve the radiative trans- 
fer equation for each bin along the vertical and 
eight inclined rays. We assume local thermo- 
dynamic equilibrium (LTE) and treat scattering 
as true absorption. 



The predicted temperature and density 
stratification as well as the dynamics of the gas 
flows are qualitatively very similar to the ones 
reported for previous simulations of convec- 
tion a t the surface of red gi ants (|Collet et al" 



a. 
id 



2 [A/B]=log(« A /«B) - log(«A/«B) G , where n A 
and fi B are the number densities of elements A and 
B, respectively, and the circled dot refers to the Sun. 



2007) and solar-type star s ( Stein & Nordlum 
1998[ lAsplund & Garcia Pered 1200 lh . Warm, 
isentropic gas ascends from the bottom and, 
as it approaches the optical surface, it begins 
to cool by photon losses. Near the surface, the 
opacity behaves as a strongly sensitive func- 
tion of temperature and decreases rapidly as 
the gas cools radiatively. As the gas becomes 
less opaque, photons escape more easily, low- 
ering the temperature even further and cooling 
the gas to the point that it becomes negatively 
buoyant and starts to sink downward. This pos- 
itive feedback causes a sudden drop of the 
gas temperature from -10000 K to ~5 000 K, 
which in the present metal-poor red giant sim- 
ulation takes place within a distance of typ- 
ically ~30 Mm or less from the optical sur- 
face resulting in very steep vertical tempera- 
ture gradients. The optical surface also appears 
to be very corrugated, with optical depth unity 
occurring over a large range of geometrical 
depths (-250 Mm). 

The temperature stratification as a function 
of optical depth for a typical simulation snap- 
shot is shown in figure Q] Superimposed is the 
corresponding stratification from a ID, LTE, 
plane-p arallel, hydrostatic marcs model atmo- 
spher e dGustafsson et al.lfl975t lAsplund et al.1 
ll997j) computed for the same stellar param- 
eters. The temperatures in the upper photo- 
spheric layers of the 3D simulation are signifi- 
cantly lower (by ~1000 K, on average) than in 
ID. This disparity can be explained consider- 
ing that in stationary hydrostatic and in time- 
dependent hydrodynamical models the mech- 
anisms controlling the energy balance in the 
upper photosphere are effectively different: ra- 
diative equilibrium, in ID, and radiative heat- 
ing and adia batic cooling due to g a s expansion, 
in 3D (see lAsplund et all 1 19991 ICollet et all 
120071) . 

We use the red giant surface convection 
simulation as a 3D, time-dependent, hydrody- 
namical model atmosphere to compute spec- 
tral lines for neutral and singly ionized Fe 
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Fig. 2. HD122563: Iron abundance derived from Fe i (left panel) and Fe n lines (right panel) as a function 
of lower excitation potential; black (grey) fil led symbols repre sent the results of the 3D (ID) analysis. The 
equivalent widths of the lines are taken from I Aoki et al.l d2007t) . 



as well as for a number of molecules (e.g., 
OH, CH, and CN), under the assumption of 
LTE. We derive Fe and O abundances by 
reproducing the m easured equivalen t widths 
of Fe i and Fe n (lAokiet al.ll2007h and in- 
frared ( IR) vibrational-ro tational OH lines at 
1.6 /urn (iBarbuy et al.ll2003t) . We then compare 
the resulting abundances with the ones derived 
with the corresponding ID marcs model atmo- 
sphere, assuming LTE and adopting a micro- 
turbulence of 2.0 kms~' for the ID spectral line 
formation calculations. 



3. Results and Discussion 

Figure|2]shows the Fe abundance derived from 
Fe i and Fe n spectral lines as a function 
of lower excitation potential. The systematic 
cooler upper photospheric temperature stratifi- 
cation of the 3D model atmosphere compared 
with the ID model result in significant differ- 
ences between the 3D and ID LTE ionization 
equilibria. In particular, the fraction of neutral 
Fe in those layers is, on average, higher in the 
3D model than in ID. Hence, at a given abun- 
dance, synthetic Fe i lines tend to be stronger 
in 3D than in ID, and a lower abundance is 
therefore required in 3D to match the observed 
equivalent widths of Fe i lines. Conversely, the 
Fe abundance derived from Fe n lines is larger 
in the 3D case than in the ID one. As a con- 



sequence, the difference between the Fe abun- 
dance values determined from Fe i and Fe n 
lines actually increases when going from ID 
to 3D. We also observe trends in the Fe abun- 
dance derived from Fe i lines with excitation 
potential in both ID and 3D, but with opposite 
slopes in the two cases. Such trends as well as 
the difference between the Fe abundance deter- 
minations from neutral and singly ionized iron 
lines are possibly due to the neglected depar- 
tures from LTE in the line formation calcula- 
tions and/or to shortcomings still present in the 
model atmospheres. 

Due to the highly non-linear temperature 
sensitivity of molecule formation, the differ- 
ences between the 3D and ID temperature 
stratifications also lead to radically different 
molecular equilibria in the two kinds of mod- 
els. From the analysis of CH lines at 4300 A, 
the derived 3D carbon abundance is ~0.2 dex 
lower than the ID value. The 3D-1D nitrogen 
abundance difference estimated from the anal- 
ysis of CN lines at 3880 A is even more pro- 
nounced and reaches ~1 dex. Figure [3] shows 
the 3D and ID oxygen abundances determined 
from the IR OH lines: again, the 3D abundance 
is lower by ~1 dex than the ID one. The 3D 
oxygen abundance values show no trend with 
lower excitation potential, which is reassuring 
as it suggests that the temperature gradient in 
the OH line formation region is probably well 
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Fig. 3. HD122563: Oxygen abundance derived 
from IR OH lines at 1.6 fjm as a function of lower 
excitation potential; black (grey) filled symbols in- 
dicate the results of the 3D (ID) analysi s. The equiv- 
alent w idths of the lines are taken from Barbuy et al. 
d2003h . 



reproduced by the 3D simulation. The ID oxy- 
gen abundance values show instead a decreas- 
ing, although relatively shallow, trend with ex- 
citation potential. 

We also determine the 3D and ID oxy- 
gen abundance from the forbidden [O i] line 
at 6300 A: we find that the 3D- ID abundance 
difference amounts to only -0.15 dex for this 
line. When we apply such 3D- ID abundance 
co rrections to th e oxygen abundances derived 
by iBarbuy et alJ ((2003) we find the corrected 
values from [O i] and OH lines to be dis- 
crepant. The reasons for this discrepancy need 
be investigated further. It may indicate that 
the temperature stratification predicted by the 
present simulations is still not accurate enough. 
We are currently re-running the metal-poor red 
giant simulation with an im proved and more 
accura te binning scheme by Trampedac h et alJ 
(2009): preliminary tests show that the result- 
ing vertical temperature gradient in the atmo- 
sphere tends to become shallower and the dif- 
ferences with the ID upper photospheric strat- 
ification appear smaller. We need however to 
wait for the simulation to fully relax before 
drawing conclusions about the consequences 
for spectral line formation and abundance anal- 
ysis. Also, one would need to investigate other 
sources of possible systematic errors such as 



the treatment of scattering as true absorption in 
the convection simulations, and, more impor- 
tantly, the possibility of departures from LTE 
in the calculation of excitation, ionization and 
molecular equilibria. 
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